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We formulate a general approach for studying the low frequency response of an interacting quan- 
tum dot connected to leads in the presence of oscillating gate voltages. The energy dissipated is 
characterized by the charge relaxation resistance which, under the loose assumption of Fermi liquid 
behaviour at low energy, is shown to depend only on static charge susceptibilities. The predictions 
of the scattering theory are recovered in the non-interacting limit while the effect of interactions is 
simply to replace densities of states by charge susceptibilities in formulas. In order to substantiate 
the Fermi liquid picture in the case of a quantum RC geometry, we apply a renormalization-group 
analysis and derive the low energy Hamiltonian for two specific models: the Anderson and the 
Coulomb blockade models. The Anderson model is shown, using a field theoretical approach based 
on Barnes slave-bosons, to map onto the Kondo model. We recover the well-known expression of the 
Kondo temperature for the asymmetric Anderson model and compute the charge susceptibility. The 
Barnes slave-bosons are extended to the Coulomb blockade model where the renormalization-group 
analysis can be carried out perturbatively up to zero energy. All calculations agree with the Fermi 
liquid nature of the low energy fixed point and satisfy the Friedel sum rule. 
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I. INTRODUCTION 

The ability to probe and manipulate electrons in real 
time constitutes one of the main challenges of transport 
in quantum dots. This program is spurred by techno- 
logical progress in guiding and processing high-frequency 
electronic signals. At low frequency, charge or spin can be 
transferred adiabatically with quantum pumps [THE] and 
single electron tunneling events can be measured by cou- 
pling the system to a nearby quantum point contact [7][5], 
the detection bandwidth is however restricted to the kilo- 
hertz regime. These low frequency experiments may con- 
trol and monitor single-charge transfer events but they 
are not able to capture the coherent dynamics of charge 
carriers. 

Early experiments at high frequency used a microwave 
source (typically above the GHz) to irradiate the quan- 
tum dot in the presence of a source-drain DC bias volt- 
age [HI [TO] • The energy of a photon hui may exceed the 
thermal energy fc^T and photon-assisted tunneling takes 
place. Electrons are then able to tunnel across the quan- 
tum dot by emitting or absorbing photons from the mi- 
crowave signal. Probing quantum dots at high frequency 
provides information on the quantum motion of elec- 
trons The microwave part of the noise power spec- 
trum corresponds to the typical energies (level spacing, 
charging energy, etc) of quantum dots of micrometer size. 
In the quantum regime Hu> > fcsT, the noise emitted by 
a nanoconductor device can be absorbed and measured 
by an on-chip quantum detector. Different schemes have 
been developed where the quantum detector, located in 
the vicinity of the source, can be a quantum dot [121 - 
[14] . a SIS tunnel junction [151117] or a superconducting 
resonator [18]. An alternative way to measure directly 
the emitted noise from the nanoconductor is to use cryo- 



genic low noise amplifiers [T9H2T] . This non-exhaustive 
synopsis of high frequency experiments demonstrates the 
vitality of research in this field. 

A fundamental and paradigmic experiment [22] in the 
topic of high frequency transport is the quantum capaci- 
tor, or quantum RC circuit. It consists of a quantum dot 
attached to a reservoir lead via a quantum point contact. 
In the experiment of Refs. [2"2~M2~3] , a single spin-polarized 
channel of the lead is connected to the quantum dot. In 
addition, the quantum dot forms a mesoscopic capacitor 
with a top metallic gate as illustrated in Fig. [T] These 
resistive and capacitive elements constitute the quantum 
analog of the classical RC circuit in series. By apply- 
ing a time-dependent voltage on the top gate, the quan- 
tum capacitor can be operated in the linear |22j or non- 
linear [2"5H2"5] regimes, and charge can be transferred be- 
tween the dot and the lead alternatively. In the linear 
regime, an AC drive changes infinitesimally and periodi- 
cally the charge of the dot. Matching the low frequency 



FIG. 1. Schematic view of the quantum RC circuit. Electrons 
coming from a metallic lead can tunnel inside a quantum dot 
where there are interactions between electrons. An oscillating 
voltage V g (t) is applied by a metallic gate coupled capacitively 
to the dot. 
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admittance of the dot with the corresponding formula for 
a classical RC circuit 

= -iujCo(l + iujCoR q ) + 0(u 3 ) (1) 

v g\ UJ ) 

allows us to define a quantum capacitance Co and a 
charge relaxation resistance R q . V g denotes the gate volt- 
age, I the current from the dot to the lead. 

Recent experiments have developed an alternative way 
of measuring the admittance of a single or double quan- 
tum dot by embedding it in a microwave cavity [251 [37] , 
see also experiments using radiofrequency reflectome- 
try [2"5H5U] . A good quality factor greatly enhances the 
coupling to the photons at the resonant frequencies of 
the cavity and the admittance can be read off the phase 
of a microwave probe. Only measurements of the ca- 
pacitance [31] were realized so far but the technique is in 
principle able to capture the charge relaxation resistance. 

The term charge relaxation resistance was coined in 
the seminal work of Buttiker, Pretre and Thomas [3"2Tf3"4"] 
where a general theory of time-dependent coherent trans- 
port was put forward. Coulomb interactions were treated 
using a discrete RPA-like model, in which the screening 
of the potential imposed by the gate is taken into account 
self-consistently. One particular outcome of this theory 
is the prediction of a quantized and universal resistance 
R q = h/2e 2 in the case of a polarized single-channel lead. 
This remarkable universality was notably confirmed ex- 
perimentally [22]. For more than one channel, R q was 
expressed in terms of statistical distribution of dwell 
times, corresponding in that case to Wigner-Smith de- 
lay times [35] [35] , and computed for chaotic or weakly 
disordered quantum dots using random matrix theory 
techniques [57^40] , Interestingly, the RPA-like screening 
approximation emerges as the leading order contribution 
in an 1 [N expansion (N being the number of channels in 
the lead connected to the dot). This expansion was de- 
vised [H] as a general method to describe the interplay 
of coherent transport and interaction in quantum dots 
driven out of equilibrium. Other aspects of the linear 
AC response of the quantum RC circuit were theoreti- 
cally addressed including the charge relaxation resistance 
and the inductive response |42j within a Luttinger model 
for a long tube connected between electrodes [43j [44] , the 
quantum to classical transition in the presence of finite 
temperature or dephasing probes [35] HB] and the effect 
of possibly strong interaction within an Hartree-Fock ap- 
proach [47H49] or by developing a real-time diagrammatic 
expansion in the tunnel coupling [50] . It has also been 
suggested that the quantum RC circuit could be used 
to detect efficiently the state of a nearby double-dot sys- 
tem [51] or to probe charge fractionalization in a quantum 
spin-Hall insulator [52] . 

The charge relaxation resistance has also been investi- 
gated for small metallic islands where the tunnel junction 
to the reservoir is described by a large number of weakly 
transmitting channels [531 154] . In this regime, a map- 
ping to the problem of a single particle on a ring subject 



to dissipation has been exploited [55] to demonstrate a 
new fixed point at large transparency associated to the 
quantized resistance R q — h/e 2 . 

Before pursuing our discussion on the linear AC re- 
sponse, let us mention that the experiment [22] on the 
quantum RC circuit was also driven into a non-linear 
regime [23H25] . A square-shaped excitation with an am- 
plitude comparable to the level spacing turns the meso- 
scopic capacitor into a single electron source. An electron 
(and a hole) is thus sent into the lead at each period of 
the driving signal. This experiment has initiated an in- 
tense theoretical activity on dynamics and quenches in 
coherent and interacting nanoscaled systems [551162] . 

The prediction of the universal resistance R q — h/2e 2 
for a single channel was recently reconsidered [631 164] by 
treating Coulomb interaction in an exact manner. In this 
way, the strong Coulomb blockade regime could be ad- 
dressed analytically, yet at the price of treating the cou- 
pling to the lead perturbatively, for either small or large 
transparency. All analytical calculations: bosonization of 
the fermionic degrees of freedom, perturbation in the dot- 
lead coupling and mapping to the Kondo model [651 166] , 
point to the fact that the resistance R q = h/2e 2 sur- 
vives arbitrarily strong interactions for all transmissions, 
thereby reinforcing its universality. In addition to this 
result, it was shown [53] that a large dot, with an effec- 
tively vanishing level spacing, also supports a universal, 
albeit different, charge relaxation resistance R q = h/e 2 . 
Interactions in the lead, for the edge state of a fractional 
quantum Hall state for instance, are added [53J [53] for 
free in bosonization. They simply renormalize the charge 
relaxation resistance R q = h/2ve 2 where v denotes the 
electron filling factor in the bulk. For v < 1/2, a quantum 
phase transition occurs [641 as a function of the trans- 
mission (dot-lead coupling) into an incoherent regime 
where resistance quantization is lost (see also Ref. [67] 
where a similar transition was obtained). Complemen- 
tary to these analytical findings, Monte-Carlo calcula- 
tions [64] have confirmed the universality of the charge 
relaxation resistance for v > 1/2 (including the case of 
non-interacting leads v = 1) , i.e. for all interactions 
and transmissions, and the quantum phase transition for 

2/<l/2. 

Surprises came from the single-channel case with spin- 
full electrons (in contrast to the fully polarized edge 
states). When a single level on the dot participates 
to electronic transport, the quantum RC circuit is de- 
scribed by the Anderson model [68] . At zero magnetic 
field and small excitation frequency, the Korringa-Shiba 
relation [69] on the dynamical charge susceptibility im- 
plies a quantized universal resistance R q = h/4e 2 , again 
in agreement with the original RPA approach [17] . Note 
that this corresponds to a weak charge response and a 
weak low frequency dissipation cx Cgi? g when the charge 
on the dot is quenched by strong Coulomb interaction 
(for example, in the Kondo regime). Finite magnetic 
fields or higher frequencies do not alter the freezing of 
charge fluctuations, but they allow processes which re- 
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distribute the spin populations on the dot and cause an 
increase of energy dissipation. The result is a giant peak 
that the charge relaxation resistance develops with either 
frequency or magnetic field [70]. Note that although the 
peak in magnetic field emerges at the Kondo energy scale 
and therefore originates from strong correlations, it does 
not contradict the Fermi liquid nature of the model at 
low frequency but arbitrary magnetic field. A Fermi liq- 
uid description [71] is thus able to reproduce analytically 
the properties of the peak, showing that the peak disap- 
pears at the particle-hole symmetric point. A generalized 
Korringa-Shiba relation can be derived that expresses the 
resistance R q in terms of static susceptibilities. As we 
shall discuss in this paper, the Fermi-liquid approach in- 
troduced in Refs. [TTJ [72] is in fact quite general and 
should apply to a variety of models. 

To summarize, analytical and numerical calcula- 
tions [63] IM] have proven that the universal resistance 
R q = h/2e 2 remains valid even for strong Coulomb in- 
teraction on the dot. Nevertheless, the physical reason 
for this universality is somehow hidden in the formal- 
ism, especially in the bosonization approach. The intent 
of this paper is to bridge the gap between the weakly 
interacting model of Refs. 32-3H HH] and the strongly 
interacting approach of Refs. |63[l64j. by proposing a gen- 
eral Fermi liquid framework that captures all interaction 
regimes within a single model, and recovers R q = h/2e 2 
in the single-channel case. 

The paper is organized essentially along two directions. 
The first part, corresponding to Sec.[TTJ presents the main 
ideas underlying the Fermi liquid approach. Due to the 
lack of phase-space available for inelastic scattering (this 
is the standard Landau argument for Fermi liquids [73T 
[75]). elastic processes dominate the physics of interacting 
quantum dots at low energy [76| . Hence, low energies are 
described by non-interacting electrons backscattered by 
the dot and, using the Friedel sum rule [77], one arrives 
at the Korringa-Shiba formula for the dynamical charge 
susceptibility. While the capacitance is proportional to 
the local charge susceptibility, Co = e 2 Xc the charge re- 
laxation resistance is expressed as a combination of static 
charge susceptibilities given by Eq. ( 24 ), with R q = h/2e 2 



in the single-channel case. The discussion is extended to 
the case of a large dot in Sec II E[ where the dot itself is 
described as a Fermi liquid of non-interacting electrons 
separated from the lead, and the charge relaxation resis- 
tance is in units of R q = h/e 2 . Both resistances, h/2e 2 
and h/e 2 , are direct consequences of the Fermi liquid 
structure Eq. ( 14 ) and of the Friedel sum rule, and are 



therefore almost independent of the gate voltage. 



The second part of the paper is detailed in Sec. Ill The 
validity of the Fermi liquid approach is explored using 
a renormalization-group (RG) analysis for two specific 
models relevant to describe the quantum RC circuit: the 
Anderson model and the Coulomb blockade model. This 
is illustrated in Fig. [2] The perturbative RG approach 
allows us, for both models, to calculate explicitly the low 
energy effective Hamiltonian in agreement with Sec. In] 



In particular, the Friedel sum rule is checked by compar- 
ing our predictions with existing results from the litera- 
ture: Bethe-ansatz calculations for the Anderson model 
and a perturbative calculation for the Coulomb blockade 
model. In addition for the Anderson model, we provide 
a rigorous mapping to the Kondo model and derive an 
analytic formula for the charge susceptibility in powers 
of the hybridization out of the particle-hole symmetric 
point. 



II. FERMI LIQUID APPROACH 

A. Hamiltonians 

The Fermi liquid approach, to be discussed below in 
the core of this section, does not rely on a specific 
Hamiltonian but rather defines a universality class for 
the low energy behaviour of coupled dot and lead sys- 
tems. Despite of this, we will introduce two particular 
Hamiltonians: the Anderson model [68] and a second 
model describing an interacting quantum dot with inter- 
nal levels that we shall call the Coulomb blockade model 
(CBM) 65, 66 . For these two models indeed, we shall 
see in Sec. |III| that they fall into the category of Fermi 
liquids at low energy. Nevertheless, one has to keep in 
mind that the conclusions of this section are in no way 
restricted to the models Eqs. ([2]), Q but would also be 
applicable for more complicated interactions on the dot 
mixing, for example, long range and short range compo- 
nents. 

The quantum RC circuit is described [70] |7T] by the 
Anderson model when the level spacing in the dot is suffi- 
ciently large and electron transport is not spin-polarized. 
The Hamiltonian takes the form 



Ham = £k<jc\ a c k< j + e d h 

a.k 

+ Ufi^hi +ty~] (cl a d a + d^cfca) j 



(2) 



with the electron operators Ck& and d a of spin a for the 
lead and the dot respectively. The lead electrons are 
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FIG. 2. Line of reasoning developed in this paper. The lin- 
ear response of interacting systems is described by a non- 
interacting Fermi liquid at low energy. It is then possible to 
compute the quantum capacitance Co and the charge relax- 
ation resistance R q . A is the cutoff energy scale in the RG 
approach. 
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characterized by the single-particle energies with the 
constant density of states v$. The electron number on the 



dot is n = 



n-j- 



with hrr — d\d„. U is the interaction 



energy, or the charging energy in the case of a quantum 
dot, related to the gate capacitance C g through U = 
e 2 /C g . Ed — —eV g is the single-level energy of the dot. 
It is tuned via the electrostatic coupling of the quantum 
dot to the metallic gate, t is the amplitude for electron 
tunneling between the dot and the lead. We will later 
need the hybridization constant T = irvot 2 . 

The second model, which we call the Coulomb block- 
ade model (CBM) [65] [66] , is appropriate for a larger dot 
with at least a few energy levels in the dot relevant for 
transport. Its simplest version includes spinless electrons 
and a single channel in the lead but it can be straightfor- 
wardly extended to N channels. The Hamiltonian splits 
as [53 ES] H CBM = H + H c + H T with 

H a = 2_j £k c L c *w + X! £l d \<j dl ° + £d ( 3a ) 

k,a J, or 

H c = E c h 2 , H T = tY, (d\ a c ka + clji^ , (3b) 

k,l,a 

where the three terms describe respectively non- 
interacting electrons with single-particle energies £k 
(lead) and £; (dot), the charging energy due to strong 
Coulomb repulsion in the dot and the tunneling of elec- 
trons between the dot and the lead. E c — e 2 /2C g and 

(7 = 1, . . . , N. Again h = J2 a n a = YjH,o d \o d ^ is the 
total number of electrons in the dot and Ed = —eV g is set 
by the gate voltage. The level spacing in the dot is finite 
but can be sent to zero for a large enough quantum dot. 
In this paper, the density of states vq is chosen for sim- 
plicity to be the same in the lead and in the dot but none 
of our results are affected by releasing this constraint. It 
will be convenient later to use the dimensionless conduc- 
tance g = N(vot) 2 . For both models, the total number 
of electrons, that is a constant of motion, is written N t . 

Although the derivations of Sec. |III| are valid only 
away from charge degeneracy, the Anderson and the one- 
channel Coulomb blockade models are both Fermi liquid 
for all gate voltages s&. Our effective Fermi liquid ap- 
proach, which predicts the resitances h/2e 2 and h/e 2 , 
is therefore applicable for all values of Ed including the 
Coulomb peaks at the charge degeneracy points. In the 
multi-channel Coulomb blockade model, the Fermi liquid 
approach breaks down only right at charge degeneracy 
(non- Fermi liquid fixed points). 



B. Non-interacting electrons 

Before addressing the general Fermi liquid approach, 
that is the central part of this paper, it is instructive to 
shortly review the non-interacting case [321434] following 
the discussion of Ref . [H] . First the coupling to the gate 
voltage £d(t)n can be gauged out by a simple unitary 



transformation U(t) = e 1 $ dt ^d(t)N t shifting a \\ single- 
particle energies by —Ed, notably Ek — > £k — £d for the 
single-particle energies in the lead. 

Coherent electrons are described by delocalized wave- 
functions propagating throughout the quantum dot. Dif- 
ferent trajectories, corresponding to single or multiple 
reflections at the quantum point contact opening the 
dot, interfere by adding their amplitudes. In the ab- 
sence of Coulomb interaction, electrons behave very much 
like photons traversing a dispersive medium: an energy- 
dependent phase shift $(£& — £d) is accumulated after 
passing through the dot. Although the discussion here 
is meant to be general, we can illustrate with a specific 
example. The phase shift reads 

„ _ gi27re/A 



\ — re i2ire/A 



(4) 



for a quantum dot embedded in a quantum Hall edge 
state, r being the reflection coefficient of the quantum 
point contact, A the level spacing in the dot and 2we/A 
the phase accumulated after a single turn around the dot. 
With this picture in mind, the current from the dot to 
the lead can be computed from the Landauer-Buttiker 
scattering formalism [32 - 34, 49 . At low frequency, 



lit) 



[V g (t)-V g (t-r)] 



(5) 



t is the Wigner-Smith delay time, or the typical dwell 
time of an electron in the dot (sf is the Fermi energy) , 



= h 



d${£ - E d ) 
ds 



= -h 



d^{£ F - E d ) 

d£d 



(6) 



and Eq. J5j) holds as long as lot <C 1 . 

Eqs. (|5) and ^ show that a non-zero current is an 
effect of the dispersive cavity. Physically, a time varying 
gate voltage £d(t) implies that different times of elec- 
tron arrivals correspond to different energies e — Edit) 
and therefore different phase shifts. All electrons are not 
slowed down in the same way and charge can accumu- 
late either in the dot or in the lead. Fourier transform of 
Eq. ([5| is compared to Eq. and yields 

h 



Ra 



Co 



2e 2 



h 



(7a) 



(7b) 



such that t = 2R q Co. The coherent regime thus gives a 
quantized and universal resistance R q = h/2e 2 . 

The same reasoning applies to the case of N channels 
in the lead connected to the dot. The result for the charge 
relaxation resistance reads [3S1 27] 

I V^N 2 t 2 

A E g =l T a = A £*=i V A 
2e 2 f^N \ 2 2e 2 



Ra 



f^N y 2e z f^N 



(8) 



where r CT is the Wigner-Smith delay time in the channel a. 
v^c = r a /h [3S] denotes the density of states for channel 
a in the dot. 
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C. Physical picture 

Simple arguments can be given to argue that the re- 
sults Eq. (7j) extend to the general case of interacting 
electrons. Despite its apparent simplicity, the expression 
Eq. ([6| for the time delay r that controls the low fre- 
quency transport, conveys a fundamental message: only 
electrons in the immediate vicinity of the Fermi surface 
participate to AC transport at low frequency. These elec- 
trons however, following the conventional Fermi liquid ar- 
gument [73-75], are somehow protected against interac- 
tions due to the restriction of available phase-space close 
to the Fermi surface. Hence, even in the presence of 
strong Coulomb repulsion on the dot, low energy elec- 
trons in the lead behave essentially as if they were non- 
interacting and the line of reasoning detailed above gener- 
alizes to the interacting case. This generalization is possi- 
ble only in the generic situation of Fermi liquid behaviour 
at low energy while non-Fermi liquid fixed points usually 
require a delicate tuning of coupling constants [7M5T] . 
We note in passing that a similar restoration of phase- 
coherence caused by Coulomb interaction was explored 
in a quantum dot T-junction geometry [76] . 

Relatedly, Eq. (7b) for the capacitance or the static 



response of the dot, coincides with an exact expression 
as we shall now argue. In the interacting case, the 
phase shift 5(e,j) = $(e/c — £d)/2 (there is factor 2 dif- 
ference between the phase shift definition used in the 
non- interacting case for instance in Refs. [22] S§]> an d 
the conventional phase shift of the Friedel sum rule [77] ) 
is related to the occupancy of the dot via the Friedel sum 
rule S/ir = (n). Substituted in Eq. ( 7b I, one finds 



Cn = -2tt 



d(n) e 2 



(9) 



where Xc is the static charge susceptibility of the dot. 
This result coincides with the mere definition of Co if we 
recall that / = dt e(n) and consider the time- integrated 
version of Eq. (JlJ in the static limit. 

Let us now clarify an important point: the similari- 
ties between the scattering properties of non-interacting 
and interacting electrons do not imply that interactions 
have no effect as they can strongly renormalize the en- 
ergy dependence of the phase shift. For non-interacting 
electrons, a comparison between Eq. (7b) and Eq. ^ 



suggests that r/h = i>do = Xc, clearly a sensible result. 
Indeed, for free electrons, the density of states repre- 
sents the number of one-particle states that fall below the 
Fermi surface upon a unit increase of the Fermi energy or, 
equivalently, a unit decrease of the reference energy Ed- 
The Pauli principle then implies that it is also the num- 
ber of added electrons, that is the charge susceptibility. 
However, in the presence of interactions between elec- 
trons, Vdo 7^ Xc in general. Xc is sensitive only to charge 
excitations while all excitations contribute to the density 
of states Vda- This difference is exemplified [70l [71] by 
the Anderson model in the Kondo regime where the lo- 
cal density of states exhibits a peak at the Fermi energy, 



mostly due to spin-flip excitations, whereas the charge 
susceptibility is suppressed by Coulomb blockade [82] . 



D. Effective model and Korringa-Shiba formulas 

We begin by some general remarks on dissipation for 
the quantum RC circuit. Integrating Eq. ([I]) with respect 
to time, we find the low frequency expansion [63] 171] 



e 2 (nM) 



C + iuClR q +O(uj 2 ), 



(10) 



for the charge on the dot in the presence of a time- 
dependent gate voltage. The comparison with standard 
linear response theory allows for deriving the correspon- 
dences 



Xc(lo = 0) 



Co 



Imx c (w) 



with the dynamical charge susceptibility 

Xc(t-t') = ^9(t-t')([n(t),n(t')})- (12) 

Following standard linear response theory, the power dis- 
sipated in the presence of an AC drive of the gate voltage, 
£d(t) = + £u cosut with e w small enough to be in the 
linear regime, is given quite generally, and specifically in 
the models Eqs. Q and by 

V=leluIm X c(u)- (13) 



The only requirement for Eq. ( 13 ) to hold is that Ed ap- 



pears in the Hamiltonian only through the term Ed n. 

Having exposed in Sec. |II C| the physical reasons for 
which the universal resistance is insensitive to arbitrary 
interaction in the dot, we proceed and develop a low en- 



ergy effective model that shall prove Eq. (7a) explicitly. 



General arguments are sufficient to reconstruct the struc- 
ture of the low energy effective Hamiltonian. Neverthe- 
less one has to keep in mind that the whole discussion 
that follows is based on the assumption of a Fermi liquid 
infrared (IR) fixed point. 

The first piece of the Hamiltonian is a free part H = 
J2k £ kc\cu- It describes the bulk lead electrons. The dot 
having a finite spatial extension, it is not able to alter 
the bulk properties (like the Fermi velocity or the effec- 
tive mass) of the lead electrons. Operators perturbing the 
free term Hq can be classified, in the RG sense, according 
to their relevance. Again the finite size of the dot im- 
plies that these operators involve only the field operators 
%jj(r = 0) = J2k c fc' r = being the entrance of the dot 
or the position of the quantum point contact, tfj ' (0) and 
derivatives thereof. There is only a single marginal op- 
erator, all other operators being irrelevant. Keeping the 
former and discarding the latter, the low energy Hamil- 
tonian assumes the form 



■ k c k >. 



(14) 



k,k' 



G 



The second term in this Hamiltonian describes a struc- 
tureless scattering potential (oc 5(r)) placed at the dot- 
lead boundary or entrance of the dot. For the two initial 
models Eqs. ^ and ([3]), the Friedel sum rule relates the 
phase shift of this scattering potential to the mean occu- 
pation of the dot [77] , namely 



■ arctan [i^v^K^d)] 



(15) 



Drawing on the ideas of Refs. [72 [72], the strategy that 
we shall adopt to compute the charge relaxation resis- 
tance is the following: the power dissipated by the gate 
voltage Ed(t ) = e d + e u cos cot in the linear regime, given 
by Eq. ( |13[ ), can also be computed from the low energy 
model Eq. (fl4|. The identification between the general 



expression Eq. ( 13 ) and the low energy calculation allows 



us to derive a Korringa-Shiba formula for Imx c (w) which 



further determines R q . 



Let us now start from the low energy Hamiltonian and 
expand it to first order with respect to (Kq = K(e d )), 



fe.fc' 



+ K'(e d ) Euj cosujt 2J c{c k ' 



(16) 



k,k' 



The scattering by the potential K 5(r) is a single-particle 
problem which can be readily diagonalized. Rewriting 
the Hamiltonian in terms of the corresponding scatter- 
ing states, characterized by the fermionic operators c k , 
absorbs the second term into the first one in Eq. ([16 1. 



Additionally, this change of basis introduces a multiplica- 
tive constant [53] in the third term of Eq. ( 16 ) and 



H = ^e k c\c k 

k 



K'(e d ) cos ojt 
1 + (irv K y 



J2 d l dk ' 

k,k' 



(17) 



This last expression is conveniently written in terms of 
the static susceptibility Xc = —d(n)/ded, which is ob- 
tained through the derivative of Eq. ([15]) 



Xc 



voK'(e° d ) 
1 + (tt^o) 2 

such that Eq. ( fl~7[ ) becomes 

X 

j£kc k c k -t- 

k 



H = e k c\c k + — cos ut V* ctcfe' 



(18) 



(19) 



in agreement with the above definition of the static 
charge susceptibility. 

Now that we have derived a more compact low energy 
Hamiltonian, we rely again on linear response theory in 
order to compute the power dissipated upon exciting the 
gate voltage sinusoidally. It involves the operators A — 
{Xc/vo)J2 k h' ^k^k' coupled to the AC drive in Eq. (19) 



(21) 



where xia(*- = iO(t ~ t')({A a (t), A a (?)}). Per- 
haps unsurprisingly, energy dissipation occurs at low en- 
ergy through the production of single electron-hole ex- 
citation. ImXA(w) is easily computed at zero tempera- 
ture for non-interacting electrons (first term in Eq. (19)), 
ImX^(u;) = hirXc^, similar to the result of a Fermi golden 
rule calculation. Comparing Eq. (H3j) and Eq. (21), wc 
obtain the Korringa-Shiba formula 



(22) 



This result, substituted in Eq. ( 11 ), recovers the universal 
charge relaxation resistance R q = h/2e 2 . 

So far, our analysis has concentrated on the single 
channel case but its generalization to N channels is 
straightforward. One then finds a generalized Korringa- 
Shiba formula 



ImXcM = hiruj ^2x1 



(23) 



The static susceptibilities Xca = —d(n a ) /ded measure 
the sensitivity of the occupations of the dot, for each 
channel a, to a change in the gate voltage. The charge 
relaxation resistance thus takes the form 



h T N y 2 

U « - 2e 2 ( N 

2jct=1 Xc 



(24) 



and resembles very much the non-interacting one Eq. (Js|) . 
They coincide after the identification v§ a = Xca, already 



discussed in Sec II C and valid only for free electrons. 
Therefore Eq. (24) gives the correct generalization of 
Eq. ([8| to the interacting case. 



E. Effective model for a large dot 



It can be checked for consistency that this low energy 
model satisfies the Friedel sum rule in the static limit 
u) — > 0. The scattering potential in Eq. ( 19 1 adds the 
phase shift 5i = —hvq^Eu to lead electron wavefunc- 
tions. The Friedel sum rule then translates this phase 
shift into a shift in the occupation number on the dot 



5(h) = 5t/ir = -Xc 



(20) 



The case of a large quantum dot deserves a specific 
discussion. By a large dot, we mean that the single- 
particle spectrum can be treated as continuous in the 
dot such that energy dissipation takes place both in the 
lead and in the dot. 

In fact, finite dots have a finite level spacing A and 
one may wonder at which energy scale the spectrum can 
be considered as continuous. One solution, proposed in 
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Ref . [53] , is to send an AC signal with a bandwidth larger 
than A in order to smear the discreteness of the spec- 
trum. In that case, the frequency uj has to be larger 
than A. A second possibility is to use a frequency w 
larger than the energy a/A^th, where -Eth > A is the 
Thoulcss energy, or inverse time of diffusion through the 
dot. Above this energy, it has been shown [S3] that the 
one-particle density of states loses its discreteness due to 
electron-electron interactions. 

The Fermi liquid picture still applies to the large dot 
with the subtlety that dot and lead constitute two sep- 
arate Fermi liquids. At low energy, the transfer of elec- 
trons between the dot and the lead is energetically pro- 
hibited [74] and electrons are effectively fully backseat - 
tered at the boundary between the dot and the lead. 
The effective model for the lead is the same as above 



(Eq. (14)) with the potential scattering strength related 
to the mean occupation of the dot via Eq. ( 15 ). 



At first glance, it may seem that the dot simply adds an 
additional channel for dissipation so that Eq. ( 24 ) would 
apply with N = 2. However, as we shall see below, there 
exists a lead/dot symmetry which reestablishes univer- 
sality in the charge relaxation resistance. We first note 
that the fact that the charging energy is ascribed to the 
dot is physically sensible but mathematically arbitrary. 
Using the fact that the total number of electrons in the 
system N t is conserved by the Hamiltonian, one can re- 
place h — Nt — hi,, hi being the number of electrons 
in the lead, and transfer the Coulomb interaction to the 
lead. Therefore, the low energy model for the dot is the 



same as for the lead, namely Eq. ( 14 ) , but the strength of 



the scattering potential for dot electrons, noted K(i ot (e d ) 
is now given by 



N t - (n) = (h L ) = —- arctan Wv Q K Aot (s d )] . (25) 
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and the charge relaxation resistance R q — h/e 2 . The 
extension to N channels is straightforward. 

III. RENORMALIZATION 

A. Outline 



The aim of this section is to justify in more detail the 
low energy form Eq. ( 14 ) for the Anderson Eq. ([2| and 
CBM Eq. ^ models. We recall that the discussion of 
Sec. II D rests on two fundamental assumptions: (i) an 



infrared (IR) Fermi liquid fixed point and (ii) the Friedel 
sum rule Eq. (15 1 applies. The two models will be dis- 



cussed separately. 

The Friedel sum rule and the Fermi liquid proper- 
ties at low energy are well-established for the Ander- 
son model [77l [85] . Our motivation is thus a practical 
one: we calculate the scattering potential K(e d ) per- 
turbatively to second order in T in the Kondo regime. 
To this end, a RG treatment is first carried out on the 
Hamiltonian Eq. ([2]) and stopped at intermediate ener- 
gies Tk <C A -c U, where A denotes the running energy 
scale (or cutoff). The RG approach remains perturba- 
tive up to these energies and can be performed explic- 
itly. At this stage, charge excitations have been com- 
pletely integrated out and the Kondo model is obtained, 
with exchange and scattering potential terms. Proceed- 
ing towards lower energies, the RG procedure becomes 
non-perturbative across the Kondo temperature Tk and 
we thus rely on the work of Cragg and Llyod [55H55] , 
They showed that, whereas the exchange term disap- 
pears at low energy, leading simply to a ir/2 scattering 
phase shift, the scattering potential is unaffected by the 
Kondo crossover up to small corrections that are negligi- 
ble within our second order calculation. The main steps 
of the RG procedure up to the IR fixed point are summa- 
rized in Fig. [3j We finally obtain the scattering potential 



An alternative formulation of the same physics is that 
the phase shift accumulated after backscattering at the 
boundary is S(ed) for lead electrons and St — 5{e d ) for dot 
electrons, where St — N/tt. Following the same steps as 
in Sec. |II D[ one finds the effective low energy Hamilto- 
nian 



Integration of the 
high energy charge modes 



Anderson 
model 



Kondo model 

including 
pot. scattering 




Fermi liquid 
including only 
pot. scattering 







k,a=L/D 



Xc 



£ w COS^i ^2 (fklCk'L ~ Ad^'d) 
k,k< 



(26) 



where L/D stands for lead/dot electrons. The Korringa- 
Shiba formula is then 



(27) 



FIG. 3. Renormalization-group (RG) analysis of the Ander- 
son model. 



v K(e d ) = - 



r 2e d + U 
2tt e d {e d + U) 

TU 2 



1 - 



TU 



TTE d (e d 



TTe d (e d + U)(2s d + U) 



hi 



U) 

u 



(28) 



This expression can be used to compute the static charge 
susceptibility \c at zero temperature, 



(Ed + U) 2 



2T 

n 



1 



1 



s d + U 



— I +2 



1 



1 



In 



-£<j 



(29) 



in agreement with a Bethe-ansatz calculation at the 
particle-hole symmetric point |89) . Let us finally men- 
tion that the preliminary perturbative renormalization 
(stopped at A ^ T K ) recovers Haldane's formula [551 |9"U] 
for the Kondo temperature of the Anderson model, 
namely 



2tt 



a -ns d (e d + U)/2UT 



(30) 



where C = 0.5772 is the Euler constant. e 1/4+c /(27r) = 
0.364 is in agreement with Ref. |90j . 

The situation is different for the Coulomb blockade 
model which is known to be a Fermi liquid at low energy 
except close to the charge degeneracy for N =^ 1 [9"T] . 
For this model, perturbative renormalization can be per- 
formed explicitly down to low energy. In this paper, this 
is done to second order in the conductance g, i.e. to 
fourth order in the tunneling matrix element t. We are 
thus able to check the scenario described in Sec. Ill El of 
two separated Fermi liquids at low energy with scatter- 
ing phase shifts in agreement with the Friedel sum rule. 
To leading order in a large N calculation, we find for the 
backscattering of lead electrons at the boundary 



K{e d ) = v t 2 In 



gc + £d 



{A[e d ]-A[-e d ]), (31) 



with A[ed] reported in ( |90| ). Translated into the dot oc- 
cupancy with the help of the Friedel sum rule, this result 
coincides with a direct calculation by Grabert [75J . 

Below we detail the perturbative calculations of the 
scattering potential K{Ed) leading to Eq. (28) and 
Eq. (31 1. A refined quantum field theory approach is not 



necessary for the leading order (first order in T) and we 
simply use the unitary Schrieffer- Wolff transformation to 
compute K(ed). This is done in Sec. |III B[ In Sec. |III C 



and |III D| the next order is obtained within a more gen- 
eral field theoretical approach to renormalization. 



B. Schrieffer- Wolff unitary transformation 

When the charging energy, U or E c , largely exceeds 
the hybridization to the lead due to tunneling, the dif- 
ferent charge states on the dot become well-separated in 
energy. For temperatures T <C U,E C , one charge state 
defines the low energy sector, the others being only virtu- 
ally occupied. The Schrieffer- Wolff unitary transforma- 
tion [93] accounts for these virtual states by integrating 
them perturbatively into an effective Hamiltonian acting 



in the low energy sector. The Schrieffer- Wolff transfor- 
mation was initially devised for the Anderson model, it 
shall be applied in this paper also to the Coulomb block- 
ade model Eq. 

1. Anderson model 

We focus on the Kondo regime defined by < — Ed < 
U, that is in-between the Coulomb peaks, with |ed|/r 3> 
ln(/7/r) in order to neglect the renormalization [55] in 
the position of these peaks. The low energy sector then 
corresponds to a single electron on the dot. The tunnel- 
ing term in the Anderson model pj, that we call Ht, 
couples subsequent charge states and is linear in t. The 
idea of the Schrieffer- Wolff unitary transformation e lS is 
to cancel this linear tunnel coupling, thereby producing 
couplings between charge states that have higher orders 
in t. This strategy is realized with the choice 



[S, ^am] 



(32) 



where H < ^ M is Eq Q with t = 0. The rotated Hamilto- 
nian assumes the form 



^AM 



JS 



H 



-iS 



ztO 

-"AM 



(33) 



After this transformation and specifically for the An- 
derson model in the Kondo regime, the coupling of the 
single-charge sector to other charge states starts as t 3 at 
most. Hence, for the purpose of a calculation up to sec- 
ond order in t, the single-charge sector decouples from 
other charge states. The resulting effective Hamiltonian 
is the Kondo model which includes a potential scattering 
term 



H'am = H a + JS ■ s 



kk'tr 



(34) 



where S denotes the spin operator for the single electron 



E 



„t 



kk' aa' ka 2 



-Ck'a-' t ne local spin of 



in the dot and s 

lead electrons. t ct(7 < is the vector composed of the Pauli 
matrices. The coupling constants are given by 



v Q J = v Q J 

V K = VqKq 



2T ( 1 

"V \e d + U 

r / 1 



i 



+ 



i 



2ir \e d + U e d 



(35a) 
(35b) 



Notice that K = at the particle-hole symmetric point 
e d = —U/2, leaving exclusively the Kondo interaction 
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with J = 8i /[/. In deriving the Kondo model Eq. (34 1, 



one discards the energy dependence of the coupling con- 
stants Eqs. (35) because the model Eq. (34) is applicable 



only for energies much smaller than the charging energy 
£ -C U, The Schrieffer- Wolff transformation is thus 
an economical way of renormalizing the Anderson model 
up to the energy scale A with Tk <C A <C U, |e<j|. 

Controlling the low energy behavior, the IR fixed point 
(A — > 0) of the Kondo model with potential scattering 
Eq. (34), has been identified by Cragg and Llyod us- 



ing a combination of analytical and numerical calcula- 
tions [8"6H8"B] . It is a Fermi liquid in which lead electrons 
at the Fermi level acquire a phase shift with three contri- 
butions: 7r/2, the phase shift corresponding to K(sd) and 
a correction proportional to T 3 that we can legitimately 
neglect in our first order calculation and even in the sec- 
Absorbing the 7r/2 



IIIC 



ond order calculation of Sec. 
phase shift into a redefinition of the lead electrons, we 
indeed obtain the low energy form Eq. (14) where K(sd) 
is given by Eq. ( 35b ) . 



With the help of the Friedel sum rule Eq. (15), the 



dot occupancy is computed from Eq. ( 35b I and then the 



static charge susceptibility at the particle-hole symmetric 
point reads 



Xc = 



d(h) 



e d = -U/2 



sr 



(36) 



in agreement with a Bethe ansatz calculation expanded 
to leading order in T/U [89"] . 



2. Coulomb blockade model 

In the Coulomb blockade model of Eq. ([3]), the tun- 
nel coupling between the subsequent charge states is also 
linear in t and the principle of the Schrieffer- Wolff trans- 
formation remains the same as in Sec. IIII B II with the 
choice Eq. ( 32 ) and the rotated Hamiltonian Eq. ( 33 ) 



where the subscript AM is replaced by CBM. 

Following Grabert [75JI22], we define the operator h, 
which gives the number of electrons in the dot, as being 
independent from the fermionic degrees of freedom di a 
and we note |n) the charge state with n electrons. In 
this representation, the tunneling term in Eq. f3]) reads 



H T =t 



n.k.l 



d\c k \n + 1) (n\ + c{di \n - 1) (n|) . (37) 



For the sake of simplicity, but with no loss of general- 
ity, we restrict ourselves here to the single-channel case 
N = 1. The operator S that generates the unitary trans- 



formation, solution of Eq. 
where 



(32), takes the form 



-it s kln c\di \n 



k,L; 



1 



Skin — 



ei -e k + E c (2n- 1) +e d 



(38a) 



(38b) 



In contrast with the Anderson model, a t 2 coupling be- 
tween each state |n) with \n ± 2) is produced by the 
Schrieffer- Wolff transformation. Nonetheless, it is possi- 
ble to remove this coupling by applying a second unitary 
transformation whose specific form will not be needed 
here. After that, the Hamiltonian becomes block di- 
agonal in the charge states up to second order in t. 
For —E c < Ed < E c , the charge state n = defines 
the low energy sector and the rotated Hamiltonian is 

ffcBM = Hq + -Hb, 

t 2 / 

ft= 2 X! [skmd\,c k ,c\di - s k ixc\did\,Ck< + h.c. 

kk'W 

(39) 

The normal ordering of operators is necessary to classify 
them according to their RG relevance. For example, 

d]c k c\,d v = 5w9(-£ l )c k c\ l +8 kk ,Q{e k )d\d l >+ : d\c k c\,di, : 

where 9(e) is the Heaviside function and : . . . : denotes 
normal ordering with respect to the Fermi sea. The last 
term in this expression describes interaction between lead 
and dot electrons; it is irrelevant at low energy and can 
be discarded. The two quadratic terms give marginal 
terms in the Hamiltonian corresponding to backscatter- 
ing of electrons at the lead-dot boundary. If we take for 
instance the first one, its contribution to Hb is 



i 2 



} J Cfe'4 ^2 6n>9(- £ i) s kio 
de, 



kk 



w 

2 r 



kk' 



i — e k — E c + Ed 

E k + E c -£ d + Dp 
Ek + E c — Ed 



(40) 



Since we consider the low energy properties of the model, 
it is consistent to take Ek = in this expression, the dif- 
ference being again irrelevant in the RG sense. Collecting 
all marginal terms in Eq. (39) and performing the inte- 
grals, we obtain the effective low energy model to leading 
order in t 



H'cbm — Ho 



9_ 



In 



E r . 



. IV kk' 



(41) 

in which the intermediate cutoff Dp has disappeared |66j 
and the dimensionless conductance g = (vp t) 2 has been 
introduced. This expression confirms the scenario pro- 



posed in Sec. HE for a large dot in which the two Fermi 



liquids, lead and dot, completely separate at low energy 
while they experience potential scattering terms with op- 
posite amplitudes. 

The phase-shift acquired by lead electrons from 
Eq. (Ell is 



5 = irg In 



E c — £d 
E c + £ d 



(42) 
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and, applying the Friedel sum rule, one finds the dot 
mean occupancy 



(n) = - = g In 

7T 



(43) 



in agreement with previous direct perturbative compu- 
tations [53 [7S1 [52] to order g. The above result Eq. Q 
is unchanged for N channels as long as g = N(yot) 2 . 

As we shall show below, the different conclusions ob- 
tained in this section carry over to the next order in g. 
We were not able to iterate the Schrieffer- Wolff trans- 
formation to next orders. We shall instead use a field 
theoretical approach to derive the renormalization of the 
two models to second order in T and g, the Anderson 
model in Sec. [mC] and the CBM in SecHEDl 



slave-bosons in this paper differs substantially from most 
works based on slave-bosons. In fact, we are interested 
in the weak coupling regime (above the Kondo temper- 
ature) and not in describing the Kondo crossover. The 
precise identification between operators is given by 



da = b\ f a + af_ a b 2l 



(44) 



ensuring the anticommutation relation did 



•44 



0. 



Clearly this operation enlarges the Hilbert space by 
adding unphysical states. The set of physical states is 
thus recovered by imposing the constraint 



(45) 



which commutes with the Hamiltonian. In this basis, the 
Hamiltonian takes the form Ham — Hq + H c + Ht, where 



C. Field theory approach I 
Renormalization of the Anderson model 

Quantum Field Theory is an efficient tool for carrying 
out RG calculations and extracting the low-energy be- 
havior. Standard diagrammatic perturbation techniques 
are not applicable to the Anderson model when expanded 
in the tunneling constant t, because of the non-quadratic 
interaction term Uh^h^. To circumvent this problem, we 



Energy 

2t d + U 



4f 
+ + 



Original 
basis 

44 io) 



41°) 



New 
basis 

4 |o> A 

n\o) x 

&Ilok 



FIG. 4. New definition of the states on the dot in Barnes' 
representation. The energies are given for an isolated dot. 

follow a representation due to Barnes [S3] and introduce 
two slave-boson operators bo and 62, corresponding to 
the states of the dot with zero and two electrons respec- 
tively. This also necessitates to define two parafermions 
operators fa destroying electrons on the dot. The rela- 
tion between the old and the new basis of states on the 
dot is summarized in Fig. [4] where |0) A denotes the vac- 
uum state of slave-bosons and parafermions. The use of 



H = J2 e kc\ a c k a + (e d + A) flf*> ( 46a ) 

ka a 

Hc = Xblb + {2e d + U + X)blb 2 , (46b) 



Ht — 



ka 



—a^kc 



(46c) 



The chemical potential A has been introduced in order to 
eliminate non-physical states by imposing the constraint 
Eq. (45). This projection is realized by setting A — > +00 



at the end of calculations [53, 94, 95 . The key point of 
the representation of Eq. ( 46 1 is that the Hamiltonian is 



quadratic in the absence of tunneling and standard dia- 
grammatic techniques become applicable. To make more 
easily contact with the Kondo model, it is convenient to 
shift the chemical potential A —> A — Ed- 
With the free propagators 



G ka( iuJ n ) =iuJ n - £fc, 
F^tyuin) =iiL> n - A, 

FQ X (iv n ) —iv n + e d - A, 



e d -U-X, 



(47a) 
(47b) 
(47c) 
(47d) 



where iuj n = (27m + l)//3 denote always fermionic Mat- 
subara frequencies and iv n = 2im/ j3 bosonic ones, the 
action corresponding to Eq. p6[) reads 
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S = - }^ c\ a (iw n )G k l(iu n )c ka (iu n ) - ^ fU iuj n )F a 1 (iuj n )f a (iuj n ) 

iu) n ka iuj n a 

- ^2 [bl{iv n )F^ 1 {iv n )b {iv n ) + b\{iv n )F^ (iv n )b 2 {iv n ) 



(48) 



ka 



The action being quadratic in the high energy modes 6o and b 2 , their integration is straightforward and gives an 
action describing a Kondo model with frequency dependent couplings 

S' = S + i [i7S CTCT / • s TT / +fcS aa >5 TT >]cl^(iu> 1 )c k > a >(iLi 2 )ft(iu n + iu 2 )U'{iv„+iu 1 ), (49a) 



J" = -2t 2 [F 2 (iv n + iui + iu 2 ) + Fo(iv n )] , K = — [F 2 {iv n + iui + iu 2 ) - F (iv n )] 



(49b) 



where So stands for the free action of lead electrons and 
parafermions. 



Free 
propagators 



moment. The action now takes the form 



k,iui n 



G k (iuj„ 



F„{iUJn) 




k, iui\ , (j 



FIG. 5. (left) Free propagators of the theory . (right) Bare 
interaction of the theory. The frequencies inside the parenthe- 
sis are the arguments of the frequency-dependent couplings in 
(49bl. 



1. From Anderson to Kondo 



We now elaborate on the connection between the ac- 



tion (49 1, strictly equivalent to the Anderson model, and 
the Kondo model. The infrared (IR) limit of the fre- 
quency dependent couplings Eq. ( 49b ) is fixed by the 

the Fermi 



poles of the Green's functions, Eq. (47a 



momentum and (47b), i.e. iv n = A, iu>i 



at 
= ioJ 2 



0. In 



this limit, the couplings J and JC reproduce exactly the 
values J and K of Eqs. (35) obtained after the Schricffcr- 



Wolff transformation, and they control operators that are 
marginal in the RG sense. Expanding the couplings in 
the frequencies iv n , ioJ\/ 2 around the IR limit generates 
irrelevant operators in the action that we neglect for the 



S = S - 



(h 



JS(r) • s(t) + K V cl a {r)c kl(J {T)h s {T) 



kk' a 



( 5 °) 

definitions S = Yj^a' ft Is 2 t f°' and ™/ = 
(s is the local spin operator for lead electrons 
defined in Sec. IIIB1). The limit A — > +oo now en- 
forces hf = 1 and, using the spin representation due 
to Abrikosov [SS], we recognize the action Eq. ((50|) as 



with the 



corresponding exactly to the Kondo Hamiltonian given 
by Eq. (34). We thus recover the same result as the 
Schricffcr- Wolff transformation in Sec. IIHB 11 

The irrelevant frequency dependences in the couplings 
Eq. (49b) die out upon lowering the cutoff A well be- 



low U, Ed and do not change the low energy Kondo form 
Eq. (50). However, along the transient region where they 



still exist, they can weakly renormalize the value of the 
constants J and K. This generates a perturbative ex- 
pansion of J and K in powers of Y/U where the leading 
order is given by the Schrieffer- Wolff results of Eqs. ( [35| . 

In order to derive this expansion, we compute the 
renormalized vertex V R (A) = Z(A)V(A) where V(A) is 
the bare vertex and Z(A) the quasiparticle weight of the 
parafermion propagator. Both quantities, to be defined 
more precisely below, depend on the running energy scale 
A. From works on the multiplicative RG approach to the 
Kondo model, it is known that the RG flow in the Kondo 
model is fully encoded in the functional dependence of 
V R on A [9"5H9T] . The strategy is thus to compute V R (A) 
in the Anderson model, or equivalently from the action 
Eq. (49), at intermediate energy Tk *C A <C U, |£<j| m 
which case the functional dependence of on A is ex- 
pected to match the Kondo one. A comparison of this 
calculation with the known expression for V R (A) in the 
Kondo model identifies the values of J, K and of the 
cutoff D that characterize the effective Kondo Hamilto- 
nian Eq. (34) (or Eq. (50)) for the Anderson model at 
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FIG. 6. Self-energy of the parafermion propagator to first 
order F. 



intermediate energy. 



2. The parafermion propagator 

We introduce the parafermion propagator T a (r — T r ) = 
~ (T T fa( T )fa( T ')) where T T denotes time ord ering. In 
frequency space, the bare propagator Eq. (47b) is modi- 
fied by the self-energy 



iui n - A - T, a (iuj n ) 



(51) 



Following Solyom's prescription [35], the vertex V(A) is 
calculated for equal incoming and outgoing frequencies 
together with the analytical continuation 

iui — > for lead electrons, (52a) 

iuj\ — > id — A for parafermions. (52b) 

A is the running energy scale of the RG flow that plays 
the role of an IR cutoff in vertex calculations, e d is the 
renormalized single-level energy of the dot obtained from 
the pole of the parafermion propagator Eq. (51 ), or 



s d = X + ReE(e d ), 



(53) 



while the residue at the pole defines the quasiparticle 
weight Z(K). Its dependence on A can be neglected, 
Z(A) ~ Z , as long as A <C max(— e d , U + Ed)- This last 
assumption is nevertheless only valid within the precision 
of our perturbative calculation (second order in T) where 
no IR divergence appears in Z(A). Close to its pole, the 
parafermion propagator takes the form 



Z 



(54) 



withz = [i-a w E(e d )]- 1 . 

The self-energy £ CT («o;„) is computed to first order in T 
to be consistent with our overall second order calculation. 
The corresponding diagram is shown in Fig. [6] with the 
expression 



S(iw„) 



* 2 E 



J2 G k (itylC( 



f(e k )+b(s d + U + X) f(e k ) + b(-X + e d ) 
iLo n + e k - e d - U - A iu n - e k - A + e d 

(55) 



The limit A — > oo is taken and the result reads, after 
summation over the energies of the lead electrons, 



A - Ed - iu n 
Do 



w , \ r / s d + U + A - iu) n 
h(iu> n ) = — In h hi 

TT \ V 

(56) 

An intermediate cutoff noted Dq has been introduced 
here to remove ultraviolet (UV) divergences. However, 
as we shall see, Dq plays no role and disappears from the 
final results. Eq. ( 53 ) for id can be solved perturbatively. 



To first order in T, one finds 



. , r /. e d + U -E d 
Sd = A + - In — h In — — 



(57) 



After standard analytical continuation, the quasiparticle 
weight is also extracted from Eq. ( 56 ) 



Z = 1 - v t 



1 



1 



s d + U s d J- L 2 Jc 



(58) 



Jo is the Schrieffer- Wolff result Eq. (35a) for J 



3. The vertex 



The vertex V(A) is given by the series of irreducible 
diagrams drawn in Fig. [7j The calculation is performed 



* — (iu>A — tw, iuj, iuj) -\- 




FIG. 7. Diagrammatic series for the vertex V. 

up to second order in T, or fourth order in t. The lead- 
ing order is the bare vertex shown in Fig. [5] which, after 
the analytical continuation of Eqs. ( [52] ) , reproduces the 
Schrieffer- Wolff expression 

V(A) = Vi = S rV • s a , a J + S TT ,S a ^K , (59) 
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with the coupling constants of Eqs. (35a) and (35b I. In 
order to derive Eq. (59), we have used that id — A to 



leading order and discarded A/max(— e d , U + Ed) correc- 
tions. 

The renormalized vertex V R (A) splits as 



V 1 



Vi + V z + V a + V\ 



(60) 



where Vz and V a,b are second order terms in V. Vz 
collects the corrections to Vi in the bare vertex of Fig. [5] 
brought by the quasi-particle weight Zo of Eq. ( |58[) and 
the renormalization of the single- level energy Eq. (|57|). It 



contains exchange and potential scattering terms 



V z =S • s^o 



J2 

4 



4K, 



In 



U) 



D 2 



J2 

iSL 

2 



JoKo 



In 



-gd(grf + U) 
D 2 



(61) 



The vertex corrections V a,b are illustrated in Fig. [7] The 
first one (a) takes the form 



V a = -— F a (iuj n )Gk(iuJA +iui — w n ) [Sp T ■ s aa J iLJn ^ lLJ , iuJA+i , 



(62) 



lUJ n — IUJ , IUJ ,ID A 



-\-iui— 4~ &/3t' ^atr' ^iuj n — iuj,iui.iuJA-\-iuj — iu n \ • 



The overall minus comes from diagrammatic rules and the Einstein convention is used for spin summation. The 
calculation of this vertex correction is carried out in Appendix [X] We find V a = Vj S • s + V K with 



V a j = -4v Q t 4 
V a K = 



1 



In- 



A 



1 



In — - In — 

(e d + U) 2 "* D ' e d (e d + U)\ n D n D 

1 



1 Vln^ + ln A 



Dn 



D J e d (e d + U) 



1 



In 



A 



1 z£^_i A 



The calculation of the (b) diagram follows exactly the same steps 

X [^*/3r ' Scr' a — iuj,iuj,icjj n -{-iw— iua ^^T^acr'^iuJA — icj,zw,io; n -(-iw— ibJAl ' 

and also contains exchange and potential scattering terms 

1 



V°, = -Avot q 



V b K = -uot 4 



e d + U _ A 
e 2 "* A ' £d(£d + E/) V £>n n D, 



1 ln A 



1 



r , 1 + In — In 

(s d + U) 2 V. D 



, e d + U A 
In — — in 



Do J ' e d (E d + U)Y" D ~D Q J ' e 2 '" D a 



1 m A 



(63a) 
(63b) 



(64) 



(65a) 
(65b) 



Adding the results from Eqs. (59) 
the renormalized vertex V R (A) = 
expands as 

^0 t2 



and (65) 
VkW) 



A 



Vj = Jo - ^Jl - VoJl In - " ) , (66a) 



m is 1 v ° j is , v ° 721 ( £ d + u 
V K = «o + —J Ko + —J In 

* 8 V ~ £ d 



(66b) 



As anticipated, D has disappeared from these final ex- 
pressions and the charging energy y/—e d (sd + U) acts as 
an effective high-energy cutoff. Remarkably, the IR cut- 
off A also disappears from the expression of the potential 
scattering term Vk where the limit A — > can safely be 
taken. However, the In A dependence in Vj signals the 



onset of the Kondo singularity which develops at low en- 
ergy and restricts the validity of Eqs. (66) to the energy 



window T K < A < ^/-e d {e d + U). 

In agreement with the leading order calculation of 
Sec. |IIIB 1| Vk vanishes at the particle-hole symmetric 
point e d = —U/2. This property is in fact expected by 
symmetry to hold to all orders in T. 



4- Kondo temperature and charge susceptibility 



In order to test our predictions Eqs. (66), we compare 
them to existing results from the literature. We show 
that Eqs. (66) give access for the Anderson model to the 



Kondo temperature and the static charge susceptibility. 
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Our expression for the Kondo temperature reproduces 
a standard result due to Haldane and the charge sus- 
ceptibility agrees with a Bethe-ansatz calculation at the 
particle-hole symmetric point. These successful compar- 
isons validate our approach. 



Eqs. (68), (69) in terms of the parameters of the An- 



Let us first consider the Kondo model Eq. (34) char- 
acterized by the exchange coupling J and the high- 
energy cutoff (bandwidth) D, with no potential scat- 
tering. In the Wilsonian language, the RG flow corre- 
sponds to integrating the model continuously over high- 
energy states thereby reducing the cutoff from D to A 
and changing the exchange coupling constant from J to 
Vj (A). Hence two Kondo models are equivalent if 
the RG flow connects them even if their bare values of J, 
D are different. V^ ondo (A) can be calculated by different 
means, but in the field theory language with Abrikosov 
parafermions [95j US] , it is defined from the same renor- 
malized vertex as the one used in this paper and illus- 
trated in Fig. [7] The one-loop calculation gives [85J |9"8] 



Vj 



lo (A) = J- u J 2 In 



A 
D 



(67) 



We already gave RG arguments in Sec. |III C 1 show- 
ing that the Anderson model maps exactly onto the 
Kondo model for energies well below the charging en- 
ergy y/—ed(ed + U). The identification between Eq. (67) 
and Eq. (66aJ thus determines the coupling constants J 
and D for the Kondo model, inherited from the Anderson 
model, in terms of the original parameters t, Ed and U. 
With the arbitrary choice 



D = y/-e d (e d + U), 



(68) 



Eq. ( 67 ) is reproduced from Eq. ( 66a I if J = Jq — (i , o/2)Jq 
or 



v J 



2TU 



2T 2 U 2 



TTS d (e d + U) n 2 s 2 (e d + U) 



(69) 



The Kondo temperature Tk is the energy scale for which 
Vj (Tk) is of order one. Among the different defini- 
tions of Tk, we use the high temperature T ^> Tk ex- 
pansion jlOOj of the spin susceptibility (in proper units) 

X{ ' 4T\ HT/T K ) 2 (In T/T K ) 2 ) K ' 

A two-loop perturbative RG calculation on the Kondo 
model Eq. ( 34 ) gives the Kondo temperature [55] 



Tk = B Da/ vq J e 



1/voJ 



(71) 



depending on J and D, where B is a prefactor of or- 
der one. The precise value, B = ei +c /2tt where C = 
0.5772 ... is Euler's constant, is obtained by compar- 
ing Eq. ( 70 1 with a weak coupling (in J) calculation 
of the susceptibility J5SJ UOlj . Substituting J and D 
in the Kondo temperature Eq. ( 71 1 by their expressions 



derson model, the formula Eq. (30) derived by Haldane 
using a completely different approach [90] is recovered. 
Sofar in this discussion, we have discarded the poten- 



tial scattering term in Eq. ( 34 ) because it does not alter 



the RG flow and the results of Eqs. (|68j), (|69j) and (|30j) are 
still valid. The potential scattering term can in fact be 
absorbed [86 88J into a redefinition of the lead electron 
wavefunctions (with a minor negligible correction as dis- 
cussed in Sec. IIIB 1 ) such that the standard Kondo RG 



flow is recovered, albeit with electrons scattered with the 
phase shift 5 a —Ti/2 — — arctan^^oA^ed)) ~ —iri>oK(e d ) 
with respect to the original electrons. Here K(ed) = V K , 
see Eq. (65b). Using the Friedel sum rule J2a $<r — n 
the dot occupancy reads (n) = 1 — 2v^K{ed) leading, for 
the charge susceptibility, to the formula Eq. ( 29 1 pre- 
viously advertised. In the particle-hole symmetric case, 
Eq. (l29|) reduces to 



sr 



in agreement with a Bethe-ansatz calculation 



(72) 



With Eq. (29), we show that the charge fluctuations in 
the Anderson model can be calculated perturbatively de- 
spite the Kondo singularity that affects only the spin fluc- 
tuations. This is another indication of the spin/charge 
separation in the Anderson model |102j . This separa- 
tion no longer occurs in the presence of a magnetic field 
and the problem of charge fluctuations becomes non- 
perturbative. Nevertheless, the Fermi liquid approach 
discussed in Sec. [H] still applies to the finite magnetic 
field case [71] . 



D. Field theory approach II 
Renormalization of the Coulomb blockade model 

Unlike the Anderson model, the CBM does not exhibit 
logarithmic singularities and the derivation of the low en- 
ergy effective model can be performed using perturbation 
theory. Nevertheless, similarly to the Anderson model, 
an expansion of the CBM Hamiltonian Eq. ^ around 
the zero tunneling limit t = is unworkable because the 
unperturbed Hamiltonian is not quadratic and Wick's 
theorem does not apply. We thus extend the approach 
of Barnes by introducing one boson operator b n for each 
charge state with exactly n electrons on the dot. The 
projection onto the physical sector 



(73) 



is realized with the chemical potential A, as in Sec. |III C[ 
taken to infinity at the end of calculations. In this 
new basis, the number of electrons on the dot is n — 
n Hi^n an d the charging energy part of the Hamil- 
tonian Eq. ([3]) reads H c = E c ^2 n r?b^ n . The charge 
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states are coupled by the tunneling term 

H T = t Yl ( d l c k^ + > ( 74 ) 

kla 

where the operator A = b n _xb n removes one electron 
from the dot. 

With —E c < Ed < E c , n = defines the low energy 
sector and, in a second order expansion in the dimen- 
sionless conductance g = N(i/ t) 2 , only the charge states 
n = ±1,±2 are virtually occupied. We thus discard all 
charge states with |n| > 2 from the action. The action 
assumes the form S = Sq + S± with 



50 = Tr 

51 = Tr 



■E^^-E^r 1 ^ 



2 



(75a) 



E biK 1 ^ + 1 e ( c L^ + 4 c fe^ f 

(75b) 



where the trace stands for the typical integral over imag- 
inary time 



Tr [O] = / drO(r) 
Jo 



The free propagators are given by 



G 



-l 

'k 
D l l 

F- 1 



-d T - et, 



(76) 



(77a) 
(77b) 
(77c) 



where E n = E c n 2 + ean denote the bare energies of the 
charge states. 



The structure of the action Eq. ( 75 1 allows the straight 



forward integration of the high-energy fields b±i and 
b±2 as detailed in Appendix [B] The action then reads 
S" = Sq + S% where 



Si =Tr -&0-F0 lfe o + t 2 E [c\vdi<jb\F\d{ lcT ,Ck>cT>bo + d|' <T Cfc o -&5F_i4, o .,d{'<r'&0 

kk'll' 
(TO"' 

+t4 E 4 O -^^ 6 F l C fe''a"^'^'' i;1 24'' ( T''' C fc''^'''- F l4<T' C fc^' 5 



kk'k"k"'ll'l"l'' 
00' a" a'" 



E 



kk' k" k'" 11' l" l'' 

a a' a" a'" 



(78) 



This action is the starting point of our perturbative ex- 
pansion around t = 0. The unperturbed action for t = 
is indeed quadratic and standard diagrammatics is appli- 
cable. Switching to the frequency representation of the 
trace, the bare interaction of the theory can be drawn 
diagrammatically as in Fig. [8] In this section, we use 
the Wilsonian RG approach and introduce the running 
energy scale A in the decomposition <j> — (f> s + <pf where 
4> = bo,Ck,di. <j) s represents the slow degrees of free- 
dom with Matsubara frequencies ui n < A and $ the fast 
ones, uj n > A, that are integrated after expanding the 
action Eq. ( 78 ) around t — 0. In the diagrams shown in 
Figs. [9l [lOl |11[ |l2"l the external lines are slow modes and 
the internal lines fast modes. In the standard RG treat- 
ment, the integration over fast modes is realized contin- 
uously to follow the evolution of the coupling constants 
under RG. This is not necessary in our case since no IR 
divergence occurs in diagrams. The fast modes are there- 
fore integrated all at once and, since we are interested in 
the low energy limit A — > 0, the energy windows for the 
fast modes in fact extends over all frequencies. 



1. Low energy theory 

After averaging over the fast modes and reexponenti- 
ating the action, one obtains the low energy form 



S — So 



-bl(r)T - 1 bo(r)+ 



V E (4>Hv(t) - 4 ct (t)4' ct (t)) bl(r)bo(r) 



kk' (7 



(79) 

where JF stands for the full propagator of the slave-boson 
bo . The opposite sign of the lead and dot scattering terms 
originates from the ordering of the operators in Eq. (78 1. 



The pole of the slave-boson propagator (recall that Eq 
0) 

1 



jFo{iv n ) = - 



iv n -X- T, (iv n ) 



(80) 



defines the renormalized slave-boson energy Eq = X + 
Hq(Eq). Close to this pole, the self-energy Eq(w„) is 
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Free propagators 



Interaction of order t 2 



G*(iw„) = 
D f (sw„) = 
F (w„) = 




Interaction of order t 

Fi(i W + ™ -.!!,) 




FIG. 8. Diagrammatic representation of free propagators in 
the CBM integrated action (a). The diagrams are represented 
in the frequency domain. The bare interaction of order t 2 (b) 
and the one of order t 4 (c) are also shown. i$ = ifli — iu)\ — 
iv\+iQ.2— ico2+i^3~ iQ-z+itiJ^+Wi. The arrows pointing to the 
center of these vertices indicate the frequency dependen ce o f 
the high-energy propagators included in the interaction (78 1. 



regular and the slave-boson propagator takes the form 

Z 



iv n - E 



(81) 



Thus the IR fixed point 



with Z = [1 - d u Eo_{E )}- 
corresponds to iv = Eq. 

The coupling constant V 
rives from the interaction vertex between slave-bosons 
and lead fermions (—V for the dot fermions), illustrated 
in Fig. 10 to lowest order in g, taken at iv = Eq for the 
slave-bosons and ioj = for the fermions. A rescaling of 
the slave-boson field 6q — > \[Zq bo in Eq. ( 79 ) removes Zq 



in the action Eq. (79) de- 



from the propagator Eq. (81 ) and renormalizes the vertex 
V — » V R . The renormalized vertex V R is therefore the 
relevant object describing the scattering of electrons as in 
Sec. Ill C| for the Anderson model. After this rescaling, 



the limit A — > +oo simply enforces &o(t)&q(t) = 1 and 
the low energy Hamiltonian corresponding to the action 
Eq. (79) becomes 



H = H n 



k k' • 



c L c k' 



(82) 



and confirms the Fermi liquid picture developped in 
Sec. (HE), see for instance Eq. (14) with K{e^) = ±V^ 



for the lead/dot electrons. 

Sees. HHP 21 and MI D 31 are devoted to the evaluation 
of V R . The slave-boson propagator is first calculated, in 
order to access Zq and Eq, then the vertex V. We finally 
apply the Friedel sum rule to the Hamiltonian Eq. ( 82 1 to 



determine the mean occupancy of the dot and compare 
with a direct calculation. 



2. The slave-boson propagator 

For an overall calculation of second order in g, only the 
first order approximation of the slave-boson self-energy 
is needed. Including the diagrams shown in Fig. [9j it 
reads 



k, iui 



I. ILJ-2 




I, i(^2 



FIG. 9. First order diagrams for the slave-boson self-energy. 



t z 



E (iv n ) = -jp X] G k {ioJi)Di(ioj 2 ) 

kleriuji y 2 

(Fi(iu n + iux - iuj 2 ) + F-x(iv n + iLu 2 — itJi)) 

0(e 1 )6(e 2 ) ( 83 ) 



+ 



-N(v t) 2 J deide 2 

0(si)0(e 2 ) 
si+s 2 + E-i + A - iv n 



£i + £2 + Ei + A - iv n 



where all electrons energies have been summed over and 
the limit A — > oo taken. The self-energy in Eq. (83) 



exhibits a linear UV divergence. The theory is regular- 
ized with the cutoff function e~ s ^ D °. Fortunately, all 
UV divergences cancel out in the final expression of the 
renormalized vertex, as discussed in Appendix [C] an d the 
limit D — > +oo is eventually taken. From Eq. ( 83 ) , one 
extracts the values of Z = 1 + d LJ 'E (X) and the renor- 
malized energy Eq — A + So (A) to order g. 



3. The vertex 

The leading contributions to the vertex V are repre- 
sented in Fig. 1 1 Q| They are similar to the self-energy 
diagrams of Fig. [9] but without contraction of the lead 
electron lines. The first diagram gives 

t 2 

— -q ^ Di(iuj n )Fi(iv + ioj - iuj n ) 



-V t 



i - m 



'e + S 1 -E(A) ! 
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V 



VI 





FIG. 10. First order diagrams for the vertex V. 



= -v t 



2J Di(iuj n )Fi{w + iuj - uj n ) 



2 



de 



WOO 
e + E l -E(A)' 



(84) 



where the analytical continuations iuj — > and iv — >• 
A+So(A) have been carried out together with the A — > oo 
limit. Adding both diagrams in Fig. 10 and expanding 
to second order in g (fourth in t), we obtain 



deide2ds3 




(85) 



The number of integration variables is reduced by chang- 
ing variables, £1+62 — > £1 and integrating over £2, leading 
to 



V = v t 2 In 



£1 



+ Nv^t 4 



de\de 3 




ei,e s >0 
S 



(86) 



\s=±l 



We note in passing that the leading order result, i.e. 
the first term in Eq. ( 86 ) , coincides with our previous 



Schrieffer- Wolff calculation, see Eq. (41 1. An additional 
contribution to the renormalized vertex V R is brought by 
the slave-boson weight Zq, when the first order correction 
to Zq multiplies the first term in Eq. (86), namely 



-Nvh 4 In 



Ei 



>o 



We finally turn to the genuine second order diagrams for 
the vertex V and therefore V R . At this point, a distinc- 
tion can be operated between diagrams scaling as Nt 4 
and those scaling as t 4 . The calculation is greatly simpli- 
fied by keeping the former and discarding the latter in a 
large N calculation. In addition, the diagrams are classi- 
fied depending on whether they involve the six-leg vertex 
of Eq. ( 78 ) (of order 0) shown in Fig. [8j those diagrams 
being listed in Fig. 11 or the ten- leg vertex (of order g 2 ), 




V-2 




v, 



ii> + iui — (uu + iuj2 — iwa 




i0 — iu + iwi — 1102 + itua 

,♦•***•«., 



^i(ii> + iu — it 



F.,(»-ifl + ius) F,(fl> + 




FIG. 11. Second order (in <?) diagrams involving the six-leg 
vertex shown in Fig. [8] and scaling as Nt 4 . The diagrams 
corresponding to the contributions V3 and V4 are not shown, 
they are similar to Vi and V2 but come with an opposite 
sign and the change E n — > E- n , see also Appendix [C] The 
diagrams scaling only as t 4 have been discarded under the 
assumption of large channel number N. 



■>■■■■ i$ 



Fiiiif + iWi — iujz) 
F 2 (iv + iw + Swi 
- iw 2 - iw 3 ) 
Fiiiif + iuij — iuig) 




■ >■■■■ iu 



Fi{iu + iaj\ — !W2) 
F 2 [ii> + iuj + iwi 



FIG. 12. Same as Fig. |ll|b ut with diagrams involving the ten- 
leg vertex shown in Fig. 18] The diagrams leading to Vio,n,i2 
are not represented, they are similar to W,8,9i but come with 
an opposite sign and the change E n — > E- n , see also Ap- 
pendix [C] Notice V a = V91 + V92 = 2V 9 i. 



again in Fig. [SJ those latter diagrams being listed in Fig. 
[l2l We shall calculate explicitly only the vertex contri- 
bution Vi and quote the results of other contributions in 
Appendix [Cj 

The expression associated to the diagram Vi is given 

by 



Vi = N 



G k (-iui)Di{iu2)D m (-i 



iu 3 



F 2 (iv + iuj + iu} 3 )Fo(iD + icu + iui\ + icj 2 + ^3), 
and, after summing over the Matsubara frequencies ^1,2,3 
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in the limit A - 
V x = Nt 4 



hi i: 



> oo, one obtains 

f(-£k)f(£l)f(-£rn) 

(e m + Ei) 2 (ei - s k - e m ) 



dexd£2 



(89) 



ie uBa >a '{£2 + E 1 ) 2 (e 1 +E 2 )' 
The calculation of the other twelve diagrams illustrated 



in Figs. 11 and 12 are not particularly enlightening and 
follow the same line as the calculation of Vi . The different 
contributions are therefore summarized in Appendix [C] 
The summation over these twelve terms together with V' 
and V" is also performed in Appendix [C] where special 
attention is paid to the cancellation of the different UV 
divergences. In the limit N — > +oo, the hnal result for 
the renormalized vertex to second order in g reads 



V R = is t 2 In 

A[£d] = 2E C V 3 
8 (2E 2 C - 2E c e d 



E-i 
-£d ( 4tt 2 



Nvlt 4 (A[e d ] - A[-e d ]), (90a) 
2 E c 



•In" 



- £d 
E c — e d 



(3E c + e d )(E c -e d ) 



hi 



E, 



E r 



(2E C + e d ) 
E r 



hi 



2 E c + e d 



4E C + 2e d 



2Li 2 



3E C + e d 
4E C + 2e d 



AE c (2E c + e d ) 
(E c + e d )(3E c + e d 



hi 



4E C + 2e a 

E r 



(90b) 



This result, substituted in the low energy Hamiltonian 
Eq. ( 82 ) , gives access to the dot occupancy 



(n) = gin §— ^ - 5 2 (A[e d \ A[-e d ]) , 
E c + £ d 

by using the Friedel sum rule 

(n) = -(7V/7r)arctan(7r^ V R ) ~ -NvqV 1 



(91) 



The result of Eq. (91 1 coincides with a direct calculation 



of the dot occupancy [711 H2] using a different pertur- 
bative approach, which validates, at least perturbatively, 
the Fermi liquid description emphasized in Sec. [Til 



IV. SUMMARY AND CONCLUSIONS 

We investigated the dynamical response of a quantum 
dot attached to a lead. We studied the low frequency 
charge fluctuations on the dot that are related to the 
admittance of this quantum circuit in the linear regime. 
We argued that the system at low energy behaves as a 
local Fermi liquid where inelastic scattering events can 



be disregarded and lead electrons are simply coherently 
backscattered by the dot. The present work extends in a 
way the general analysis built by Buttiker and coworkers 
to study the dynamic admittance of mesoscopic conduc- 
tors, by including arbitrarily strong Coulomb interactions 
within the dot. We avoid the introduction of an approx- 
imate self-consistent potential on the dot. 

By computing the power dissipated by the external AC 
drive, we were able to derive a set of general formulas for 
the quantum capacitance and the charge relaxation re- 
sistance that characterize the response of the dot at low 
frequency. Remarkably, the results are essentially not so 
different from the weakly interacting picture. The funda- 
mental reason is that electrons close to the Fermi energy 
do not feel strong interactions as a result of Pauli block- 
ing, following the standard Fermi liquid argument [75] . 
Our approach is naturally not applicable for models ex- 
hibiting non- Fermi liquid physics such as the two-channel 
model close to charge degeneracy. In these models, the 
mere definition of a capacitance and a charge relaxation 
resistance is elusive as a result of unconventional scaling 
laws. 

Another important assumption in our work is that the 
Friedel sum rule is satisfied. This seems to be the case 
for the two models we investigated, but it is certainly not 
general, even for a Fermi liquid fixed point. However, the 
fact that dissipation comes from the time-dependence in 
the phase shift felt by lead electrons does not rely on the 
Friedel sum rule and should apply to more general mod- 
els such as double-dot geometries. Also we only consid- 
ered the case of zero temperature or temperatures much 
smaller than the charging energy. At higher tempera- 
tures, inelastic modes are excited and the analysis pre- 
sented here is not applicable. 

In order to put the Fermi liquid picture on firm 
grounds, we calculated, based on a renormalization-group 
analysis, the low energy theory for two representative 
models describing the quantum RC circuit: the Ander- 
son model and the Coulomb blockade model. In addition 
to providing an alternative demonstration of the mapping 
from the Anderson model to the Kondo model, we find 
that the phenomenological low energy theory proposed in 
this paper is recovered perturbatively for the two models. 
The Friedel sum rule was also checked explicitly. 

We conclude with a technical remark regarding 
Ref. jSHj where the case of a tunnel junction (infinite N) 
was considered. The definition of the charge relaxation 
resistance R q involves formally the limit of vanishing fre- 
quency u> — > 0. At finite temperature, this limit does 
not commute with an expansion in the tunneling g or T 
as emphasized in Ref. [S3] (see also Ref. |103j ). At zero 
temperature however, the two limits commute [104] and 
the perturbative expansion of Sec. |III| becomes justified. 
Moreover, at finite by small temperatures, a large fre- 
quency huj ^> ksT is sufficient to suppress the singular 
term oc 1/ui and the results of Sec. Ill are still valid. 

We thank K. Le Hur and L. Glazman for stimulating 
discussions. 
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Appendix A: Calculation of the vertex V 



Once considered 



1 3 
(S^ T • s acr ) (S T '0 ■ s a i a ) = —S T i T ■ s a i a + y^S tt i S a 



(Al) 



it is possible to determine which term will contribute to the exchange and to the potential scattering part of the 
vertex V a = S • sVj + V%, where 



V} = 4— ^2 F(iuj n )Gk(iwA + iu - iuj n ) [F 2 (iuA + iu) + F 2 (iua + iu)F (iu n — iu;)] 

P 



Vj C = —-p ^ F(iu) n )Gk(iu)A + icu - iu) n ) [F 2 (iu n - iu) + F 2 (iua + iu) + F {iuj n - iu)F 2 (iuA + iu)] 

We detail the calculation of the Matsubara sums only for the following example: 

/(A) 



(A2a) 
(A2b) 



~a y2 F {iu n )G k {iu A + iu- iu n )F^(iu n - iu) = i/ / de 

P C J -Do 



(iuiA + iu — X — e)(sd — w) 5 



fjiuA + iu-e) f(z) 

(iuiA + £d — A — e) 2 (io;A + «u» — e — A) dz (z — X)(iujA + iu — z — s) 



Z=i0J + \ — €d 



(A3) 



Notice f(iu>A + iu — e) = /(— g). The analytical continuations of Eqs. (52) can now be performed and the first 
deviation from A of id in Eq. (57) can be neglected to this order. The denominators do not depend on A anymore 
and, as required by the projection technique, the A — > oo limit can be taken. All the contributions of the poles which 
were proportional to A disappear and the only remaining integral is 



vq 



de 



vq 



o (£ + A)(e d -sy 



In^-lnA 
D D 



It is also possible to take the infinite bandwidth limit Do — > oo in this expression, but it is convenient to keep a finite 
cutoff for the moment. The same kind of calculations gives analog results for the remaining sums 



i ^2 F{iu> n )G k {iuA + iu — iu n )F (iu n - iu)F 2 (iuA + iu) 



V 



s d (£d + U) V D 



In — In — — - 

Do 



\ F ( iu) n)Gk{iuA + iuibj n )Fl(ibJA + iu) 



VQ 



(e d + uy 



In 



A 
Da 



(A4) 
(A5) 



These last expressions allow us to obtain Eqs. (63). 



r 



Appendix B: Integration of the high energy modes This action can be written as S = So + Si + S—i, with 
in the Coulomb blockade model 



The action Eq. ( 78 ) being quadratic in the high energy 



modes b± 2 , their integration is straightforward and brings 
separate terms to the effective action 



Si = Tr 



b 2 :t 2 ^2 Tr cladi^b^dl^Ck'a'h 



(Bla) 



kk'W <jo' 



b_ 2 :t 2 Tr Ul^bliF-scl^dt^b-i] . (Bib) 



kk 1 ll 1 ■ 



-b{®i% 

{ C k<T dl <? b O b l + 4a C k<rb\bo 

kla 

S-i = Tr [-bl^Z\b-i 



+t 



( C iA* b -l b + 4a C k<rblb. 

kla 



(B2a) 



(B2b) 
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and the new effective propagators 



^=F^ l -t 2 Aj^F 2 d\, a ,c k , a ,, (B3) 

hh'll'uu' 

*Z{=Flt-t 2 Y, dl C *<rF-2c\ la ,d V <r: (B4) 



kk'W < 



The integration over the b± x modes is still Gaussian and 
the effective action without high energy bosonic modes 
is finally obtained 



S=S + T r [-blF - 1 b + 

t 2 J2 c L d iA®i d L> c k'<r'bo+ 



kk'W ' aa' 



(B5) 



t2 dl C *A®-A'a' d l'<r ,b 



kk'W aa 1 



The $ operators can be expanded pert ur bat ively in t 



^l=Fl+t 2 Y FlAa^ F A-a'Ck-a'F Xl (B6a) 
kk'W aa' 

$_i = F_i+t 2 Y F_ x d\ a c k(T F_2c[, (Tl d v<T .F- X . 

(B6b) 



kk'll'aa' 



This gives the action Eq. ( 78 1 



Appendix C: Summation of all the contributions of 
V R in the Coulomb blockade model 



All the co ntrib utions corresponding to the diagrams of 
are listed below (with C = Nv^t A ) 



11 



and 



12 



Figs 

Vi = -C J de x de 2 
V 2 = C J de x de 2 
V 3 = C f de x de 2 



£l 



{e 2 +E 1 ) 2 {e 1 +e 2 y 

£l 

(e 1 +E 1 y{e 1 +s 2 y 

El 



C l de\de 2 



2C / de 1 de 2 



V 4 
V 5 

V 6 = 2t' / dc l d: 2 
V 7 
V 8 

v 9 

Vio 



e 2 + E_ 1 ) 2 (e 1 +e 2 ) 
£i 



{e x +E_ x ) 2 {e x +e 2 ) 
£l 



-C / de\de 2 
-C / de\de 2 



'e 1 +e 2 ){e 2 +E l ){e x +E_ x y 

£l 

{e x + e 2 ){e x + E x ){e 2 + E_ x y 

£i 



-2C J de x de 2 
C / de x de 2 



{£2+E l y{e l +e 2 + E 2 y 

£jl 

(s x +E 1 y(e 1 +e 2 + E 2 )' 

£i 



(e 2 + E x )(e x +e 2 + E 2 )(e x + E x ) ' 
£i 



V n =C / de x de 



Vi 



2C / de x de 2 



{e 2 + E_ x y{e x +e 2 + E_ 2 y 

£i 

(ei+ J B_ 1 ) 2 ( £l + £2+^-2)' 

El 

(ei + £U)(ei + e 2 + E- 2 )(e 2 + £U) ' 

(CI) 



where integrals run over the Ei 2 > domain, to which 
the contributions of Eq. ( 86 ) and Eq. ( 87 ) must be added. 



Whereas each term in Eqs. (CI I, (86) and (87) suffers 



from an UV divergence, the summation over all contri- 
butions is finite and does not depend on the cutoff proce- 
dure. We shall adopt a sharp cutoff at energy D in the 
following. Moreover, the calculation exhibits a particle- 
hole symmetry: for example, V3 can be viewed as the 
symmetric of Vi, they have opposite sign and E n ex- 
changed with E_ n , The result will then be necessarily of 
the form A[sd] —A[—£d], which implies that any constant 
independent of £d will be ignored during calculations. 
The dilogarithm function appears 



f 

Li 2 (z) = dt 

J z 



t 



and the following equalities will be exploited 
Li 2 (x) + Li 2 (l — x) = lna;ln(l — x), 



(C2) 



(C3a) 



Li 2 {x) + Li 2 



= — In 2 x — iirlnx, (x > 1). 

O £4 



(C3b) 
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As an intermediate step, we find (we omit the Nvffi 
factor) 



with the contribution of Eq. ( 86 1 



E-\ 

V' = [ de x de 2 - r 1 r ( -. £l 

a J (e 2 +E 1 ) {(ex+Ex 



(C4) 



+ 51 

) (ei+^-i) 



^ -ln£> - ^lnE + lnEi + — ilnE_i, 

(C5) 



Vi + v 2 



v 7 + v 8 



-21nEi lnE + In 2 D 

2 In D +ln 2 Ei + 21nEi, 
-Do Ea 
"E x 

1 



Ei (E 2 — Ei ) 



lnD + 

(E 1 E 2 lnE 1 -E%lnE 2 ), 



V s + V 6 + V' 



e d 2 



— In 2 ^ 1 

E c E_i 



2 In 



Ei 

E-i 



E 

V 9 = 2 In In A, - In 2 D a + — In 2 E x 

E c 

Ei ir 2 E 2 ir 2 E 2 (E 2 -Ei 

1 1 Ei 2 

E c 2 E c 6 E c \ E 2 

+ — In 2 E 2 - — In Ei In E 2 , 



where V£ is obtained from V' a by particle-hole symmetry. 
It can be checked explicitly that the terms depending on 
the cutoff Do in the above expressions cancel out when 
the summation over all contributions is carried out. One 
is left with 



Vi + V 2 = In Ei+2\nEi, 
1 



V 7 + V 8 = 



Ei(E 2 — Ei] 



{EiE 2 \nEi-E^\nE 2 ), 



V^lnEi + ^ilnE.!, 
t,i 



v 5 + V 6 + V" 

v 9 



Ei 



£ d 2 £<i , 2 Ei 

— — 7r — — In — |-21n- 

E c E c E_i E_i 

E ± hl 2 E _ ElTT 2 E2 7T 2 

E c E c 2 E c 6 



, E 2 ( E 2 — Ei 

H ^2 

E c 2 I E 2 



E2 

E r 



^\n 2 E 2 
2E C 

In Ei lnE 2 . 



Adding the particle-hole symmetric terms, one finally ar- 
rives at Eq. ([901). 
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